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玉米 籽粒 的 传 质 干 燥 模 拟 及 实验 分 析 
黄 凯 陈 兴 付 陈 鹏 扶 王 凤 资 涂 广 赵 宇 刘 相 东 杨 德 勇 
中 国 农业 大 学 工学 院 , 北京 100083) 
摘 要 : 为 了 解 干 燥 过 程 中 玉米 籽粒 内 部 水 分 的 扩散 过 程 ， 优 化 干燥 工艺 参数 ， 本 文采 用 模拟 研究 和 实验 研究 方法 
分 析 干 燥 过 程 中 玉米 籽粒 内 部 水 分 分 布 随时 间 的 变化 规律 。 玉 米 籽 粒 由 种 皮 、 角 质 胚 乳 、 粉 质 豚 乳 和 胚 四 组 分 组 成 ， 
四 种 组 分 水 分 扩散 系数 各 不 相同 ， 本 文 分 别 假设 玉米 籽粒 由 单 组 分 均 质 体 和 多 组 分 非 均 质 体 组 成 ， 分 别 建立 了 玉米 
籽粒 的 干燥 数学 模型 ， 利 用 COMSOL Multiphysics 模 块 模拟 研究 了 玉米 籽粒 内 部 的 水 分 变化 情况 ， 并 通过 玉米 籽粒 
的 薄 层 干燥 实验 进行 了 实验 验证 。 研 究 结果 表明 ， 建 立 的 两 个 模型 均 可 有 效 模拟 玉米 籽粒 薄 层 干燥 过 程 ; 80YC 的 模 
拟 值 与 实验 值 的 差异 较 70C 的 小 ; 玉米 籽粒 多 组 分 模型 精度 优 于 单 组 分 模型 精度 。 
关键 词 : 玉米 籽粒 ; 多 组 分 ; 单 组 分 ; 注 层 干燥 ， 质量 传递 
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in Drying a Single Corn Kernel 
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Abstract: To study the moisture diffusion of a single corn kernel during drying, and to optimize the 
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drying parameters, the moisture content distribution in the corn kernel over time was researched by the 
methods of simulation and experimental analysis. Corn kernel consists of four components of seed coat, 
horny endosperm, farinaceous endosperm and embryo, all of whose diffusion coefficient differ. Then two 
hypotheses are proposed: one is that the corn kernel consists of isotropic monocomponent; the other is that 
the corn kernel consists of multicomponent described above. Based on the hypotheses, different drying 
mathematical models of a single corn kernel were built, and solved by the software of COMSOL 
Multiphysics for mass transfer Thin layer drying experiment of the corn kernel was conducted to validate 
the models. The simulation and experimental results indicated that both of the models could explain the 
drying process of corn kernel well. The difference between simulation and experimental results under 
80C is smaller than that of 70C. And the simulation precision of multicomponent model is better than 
that of monocomponent model. 
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0 引言 


合理 的 干燥 工艺 对 于 干燥 行业 至 关 重 要 ， 就 玉 过 核磁 共振 、CT 扫描 等 技术 手段 只 能 确定 某 个 时 
米 而 言 ， 干 燥 工 艺 会 影响 玉米 开裂 紊 ， 储 藏 时 间 以 刻 茶 个 截面 处 的 水 分 分 布 们 ,水 分 在 不 同 籽 粒 结构 
及 总 体 产品 质量 中 ， 同 时 还 会 造成 能 源 浪费 。 目 前 ， 中 的 扩散 速度 不 同 ， 决 定 了 将 玉米 分 层 分 析 的 必要 
越 来 越 多 的 干燥 工艺 , 可 以 通过 模型 模拟 "再 进 性 。Takhar 考虑 了 玉米 籽粒 的 不 均匀 性 ， 研 究 了 
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行 一 定 的 实验 以 检验 模型 ， 最 后 调整 模型 以 改善 干 玉米 不 同 成 分 的 扩散 系数 门 。 干 燥 过 程 中 ， 玉 米 不 
燥 工艺 ， 达 到 提升 产品 质量 和 节能 减 排 的 目的 。 同 组 分 的 质地 不 一 样 ， 则 不 同 组 分 的 水 分 扩散 速度 

玉米 籽粒 形状 不 规则 ， 组 织 结构 较为 复杂 ， 在 也 不 一 样 。 同 时 ， 籽 粒 内 部 的 水 分 梯度 与 籽粒 内 部 
不 同 干燥 阶段 不 同 截面 处 的 水 分 分 布 各 不 相同 。 通 结构 息 上 息 相 关 ， 其 内 部 结构 可 以 看 作 种 皮 、 角 质 胚 
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乳 、 粉 质 胚乳 和 胚 四 部 分 组 成 的 组 合体 世 2， 其 扩散 
系数 各 不 相同 ， 而 这 都 与 籽粒 的 温度 和 水 分 含量 有 
关 。Chenl 根 据 玉 米 复 杂 结 构 ， 将 每 个 组 分 分 离开 
来 ， 分 别 测 得 了 它们 的 干燥 曲线 ， 通 过 运用 
COMSOL Multiphysics 对 菲 克 第 二 定律 进行 求解 ， 
得 出 了 不 同 成 分 的 扩散 系数 与 干燥 温度 和 含水 率 的 
关系 , 为 本 文 玉米 扩散 系数 的 选择 提供 了 理论 依据 。 
张 世 伟 IW "根据 热风 条 件 下 的 玉米 籽粒 内 部 传 热 
传 质 过 程 , 对 三 维 模型 进行 了 理论 分 析 与 模拟 研究 。 
通过 医用 CT 扫描 了 玉米 籽粒 ， 将 扫描 得 到 的 高 精 
度 图 片 通过 MIMICS 和 ANSYS 软件 优化 后 ， 建 立 
了 玉米 籽粒 的 三 维 几何 模型 ， 再 导入 COMSOL 
Mnultiphysics 模拟 软件 中 完成 求解 计算 。 这 种 方法 可 
以 用 于 更 好 地 描述 玉米 籽粒 内 部 的 质量 传递 过 程 和 
水 分 随 干燥 过 程 的 变化 情况 。 所 以 ， 本 文 拟 通 过 构 
造 多 组 分 与 单 组 分 玉米 籽粒 结构 的 物理 模型 ， 通 过 
COMSOL Multiphysics 模拟 与 薄 层 干燥 实验 进行 对 
比分 析 , 以 检验 玉米 籽粒 不 同 组 分 模型 的 模拟 精度 。 
1 干燥 传 质 模型 

1.1 玉米 籽粒 物理 模型 

1.1.1 多 组 分 几何 模型 

段 设 玉 米 籽粒 在 宽度 方向 上 结构 对 称 ， 因 此 只 
需 计算 玉米 籽粒 的 二 分 之 一 即 可 ， 其 物理 模型 如 图 
1 所 示 。 其 中 玉米 籽粒 的 四 部 分 体积 如 表 1 所 示 。 


DY \ 


图 1 PRO/E 构建 的 玉米 种 皮 、 角 质 胚 乳 、 粉 质 胚乳 、 胚 


Fig.1 Seed coat horny endosperm, farinaceous endosperm, 


embryo of corn kernel constructed by PRO/E 
表 1 玉米 粒 四 组 分 体积 
Table 1 Four components volume of corn kernel 
种 皮 ”角质 胚乳 ” 粉 质 胚乳 。” 胚 
体积 /mm 5.8 137.3 63.7 26.5 
体积 含量 比 2.5% 68.3% 31.7% 11.4% 
1.1.2 单 组 分 几何 模型 
假设 玉米 籽粒 在 宽度 、 厚 度 方向 上 结构 对 称 ， 
因此 只 需 计算 玉米 籽粒 的 四 分 之 一 即 可 ， 其 物理 模 
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型 如 图 2 所 示 。 取 玉米 籽粒 重心 为 坐标 原点 , x 轴 、 
2 轴 分 别 代 表 玉 米 籽粒 的 宽度 和 厚度 方向 。 


图 2 四 分 之 一 玉米 籽粒 物理 模型 


Fig. 2 A quarter of corn kernel physical model 

其 中 ， 玉 米 籽粒 的 多 组 分 几何 模型 与 玉米 籽粒 
的 单 组 分 几何 模型 总 体 尺 寸 一 样 。 
1. 2 热风 干燥 数学 模型 
1.2.1 多 组 分 玉米 籽粒 的 热风 干燥 数学 模型 

干燥 模型 假设 : 1) 干燥 初期 ， 玉 米 籽 粒 内 部 温 
度 和 水 分 分 布 均匀 ; 2) 不 考虑 干燥 过 程 中 玉米 籽粒 
的 收缩 变形 ; 3) 玉 米 籽粒 各 组 成 部 分 分 别 为 各 向 同 
性 的 均匀 体 ，4) 水 分 以 液态 形式 扩散 到 外 部 边界 ， 
在 玉米 籽粒 表面 汽化 。 

干燥 过 程 中 ， 籽 粒 内 部 的 传 质 方程 描述 为 : 

尝 - 至 [5 各 3pY (1) 

Ot ol po) WV WW) 0 

干燥 初始 条 件 为 : 

t=0 ,X=X, (2) 

边界 条 件 为 : 

-DT =h (XX) @) 


,= 0 
式 中 : 1 为 干燥 时 间 ，s; 为 玉米 干 基 水 分 含 
量 , kg/kg; h 为 对 流传 质 系 数 , m/s; Sh 为 Sherwood 
数 ; d 为 等 效 粒 径 , mm; XX 为 平衡 水 分 含量 , kg/kg; 
D; 为 玉米 不 同 组 分 的 水 分 有 效 扩散 系数 ， 寺 1, 2, 3， 
4, ms。 
1.2.2 单 组 分 玉米 籽粒 的 热风 干燥 数学 模型 
干燥 模型 假设 1) 干燥 初期 ， 玉 米 籽 粒 内 部 温 
度 和 水 分 分 布 均匀 ;2) 不 考虑 干燥 过 程 中 玉米 籽粒 
的 收缩 变形 ; 3) 玉 米 籽粒 为 各 向 同性 的 均匀 体 ;，4) 
水 分 以 液态 形式 扩散 到 外 部 边界 ， 在 玉米 籽粒 表面 
汽化 。 干 燥 初 始 条 件 、 边 界 条 件 与 1.2.1 一 致 。 
干燥 过 程 中 ， 籽 粒 内 部 的 传 质 方程 描述 为 : 


=v(DVX) (5) 
式 中 : DD 为 水 分 有 效 扩散 系数 ，m2/s。 

1.2.3 玉米 籽粒 的 热风 干燥 传 热 数 学 模型 

干燥 过 程 中 ， 闻 粒 内 部 的 传 热 方程 描述 为 : 


67 06/ 67 of,o7 of,o7 
0 
Ot Or\ Ox OV\ 0 OZ\ Oz 
干燥 初始 条 件 为 : 
t=0,T=7T (7) 
边界 条 件 为 : 
oT V ox 
4 h (7 1) plhs, C,(T 7)]3 A (8) 


式 中 : 7 为 玉米 籽粒 温度 ，K; To 为 玉米 籽粒 初 
台 温度 ， 开 ; 为 热风 温度 ，K; 上 为 玉米 籽粒 的 导 
热 系 数 ，W/mK)， 志 为 籽粒 内 部 温度 了 沿 外 法 线 7 
的 导数 ; 有 为 对 流传 热 系数 ，W/(m>.K);， 玉米 籽粒 的 
比热容 Cj=2010J/(kg.K) 5， 玉米 籽粒 密度 p=1150 
kg/mt 1，hjs 为 玉米 籽粒 的 汽化 潜 热 ，Jkg，F 为 玉 
米 籽 粒 的 体积 ，m ;4 为 玉米 籽粒 的 表面 积 ，m”。 


2 实验 材料 与 方法 


2.1 实验 材料 

玉米 ( 德 美亚 1 号 )， 无 破损 、 霉 变 粒 ， 色 泽 气 
味 正常 。 收获 时 初始 含水 率 0.38 土 0.02 (d.b.) kg/kg， 
收获 后 装 入 密封 伐 ， 储 藏 在 冰箱 中 (5 土 1 )C， 待 实 
验 前 取出 凉 至 室温 。 
2. 2 主要 实验 装置 

数字 型 洞 道 干 燥 实 验 装置 (浙江 中 控 科 教 仪器 设 
备 有 限 公 司 ): 鼓风机 (BYF7122) 额 定 功率 370W; 电 
加 热 器 额定 功率 45kW ; 干燥 室 尺 寸 为 
180mmx180mmx1250mm; 称 重 传感器 为 CZ500 型 ， 
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图 3 干燥 装置 示意 图 


1 一 鼓风机 ; 2 一 管道 ，3 一 进 风口 ; 4 一 加 热 器 ; 
5 一 干燥 室 ， 6 一 气流 均 布 器 ，7 一 称 重 传感器 ; 
8 一 物料 ;9 一 玻璃 视 镜 门 ; 10，11，12 一 蝶阀 

Fig. 3 Schematic diagram of drying equipment 


1-Fan; 2-Pipe; 3-Air inlet; 4-Heater; 5-Drying chamber; 
6-Air distributor; 7-Load cells; 8-Material; 9-Glass mirror 
doors; 10,11,12-Flygate 
2. 3. 2 计算 机 模拟 实验 
表 2 传 质 传 热 过 程 参数 


Table 2 Parameters in the heat and mass transfer 


process 
模型 参数 参数 设置 
玉米 表面 扩 散 系 数 D/ m/s 7.817x 105 0 
种 皮 扩 散 系数 D1/ m/s 3.76 x 10™ EE 


硬 质 胚乳 扩散 系数 Dy/ m/s 
软 质 胚乳 扩散 系数 D3 mAs 


4.08 x 107 区 
9.65 x 107e 和 


0~300g。 
其 他 仪器 设备 : DHG-9140A 电热 恒温 鼓 风 干 燥 
箱 ( 上 海 精 宏 实 验 设备 有 限 公 司 ),FB224 电子 天 平 (上 
海 终 宇 恒 平 科学 仪器 有 限 公 司 )，HM70 温 湿 度 仪 ( 芬 
兰 Vaisala 公司 )。 
2.3 实验 方法 
2. 3. 1 玉米 籽粒 热风 薄 层 干燥 实验 

设置 干燥 实验 条 件 : 干燥 初期 籽粒 内 部 水 分 、 
温度 分 布 均匀 ， 以 玉米 初始 含水 率 39.33%(d.b.)， 初 
始 表面 温度 25'C， 干 燥 风 温 70C 以 及 80'C( 注 层 干 
燥 实 验 发 现 ， 热 风 温度 越 高 ， 玉 米 籽粒 的 干燥 速率 
越 快 ， 但 并 非 温 度 越 高 干燥 效果 越 好 (1 由， 因而 本 文 
选取 70C 和 80C 作 为 热风 温度 )， 风 速 1.2 m/s，3 
气相 对 湿度 12 土 2% 为 初始 条 件 ， 玉 米 籽粒 干燥 至 安 
含水 率 14%。 分 别 测定 玉米 籽粒 的 干燥 曲线 。 
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胚 扩 散 系 数 Duym2/s 8.79x10 “eS 
舍 伍德 数 Sh 2.0+0.522Re"’ Sc’ 
玉米 等 效 粒 径 dimm 4.83 

常数 D。 2.5x10-” 


exp(—1.74—3.56X+ 


KH W/m:K) 


4.72x10 T+6.48X° — 


1.5x10°7T* +6.72x10° XT7) 

水 的 汽化 潜 热 hy J/kg 383.65x (647.1— 7) 
玉米 籽粒 的 汽化 潜 热 
ha/ Jkg 

模型 的 数值 求解 采用 COMSOL Multiphysics 
(4.3a) 中 的 对 流 扩 散 模块 和 传 热 模块 上 ， 根 据 传 执 
传 质数 学 模型 自 定 义 的 控制 方程 ， 对 玉米 籽粒 进行 
模拟 实验 ， 模 拟 实验 条 件 与 实验 条 件 相 同 。 在 
COMSOL Multiphysics 设置 过 程 中 ， 导 入 物理 模型 
后 ， 在 单元 尺寸 中 ， 校 准 选择 普通 物理 网 格 标准 ， 


h x(1+0.8953e ™) 


通过 网 格 无 关 性 验证 表明 ， 取 普通 物理 网 格 作 为 计 
算 网 格 满足 精度 要 求 。 
传 质 传 热 过 程 各 参数 如 表 2 所 示 。 


0 数 的 获取 
含水 率 计算 公式 为 : 
(8) 
W, 


式 中 : 有 刺 为 在 任意 干燥 上 时 刻 的 总 质量 ，g; 
Wo 为 干 物质 质量 ，g 

结果 分 析 与 讨论 

模拟 干燥 5min 后 ， 籽 粒 内 部 不 同位 置 处 温度 
差异 很 小 ， 温 度 梯度 近似 为 零 ， 此 后 仅 考 虑 水 分 梯 
度 对 水 分 扩散 的 影响 。 由 图 4 可 知 ， 热 风 温度 为 
70'C,， 干燥 12min 时 ， 单 组 分 玉米 籽粒 平均 含水 率 
模拟 值 和 实验 值 相差 最 大 ， 最 大 差 值 为 0.0389， 相 
对 误差 为 11.39%, 误差 来 源 可 能 是 由 于 实际 物料 与 
模拟 假设 参数 间 的 差异 。 模 型 中 采用 一 维 球体 湿 分 
扩散 模型 得 到 的 经 验 水 分 有 效 扩散 系数 ， 且 将 玉米 
籽粒 简化 成 均 质 体 ， 而 实际 玉米 籽粒 由 种 皮 、 胚 、 
角质 胚乳 、 粉 质 胚 乳 等 组 分 组 成 ， 且 各 部 分 扩散 系 
数 不 同 ， 对 其 简化 必然 引起 偏差 。 相 比 而 言 ， 多 组 
分 玉米 籽粒 平均 含水 率 模 拟 值 和 实验 值 具有 较 高 的 
重合 度 ， 最 大 差 值 为 0.0138， 相 对 误差 为 4.56%。 
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干燥 时 间 /min 
图 470°C 下 单 组 分 与 多 组 分 模拟 值 与 实验 值 干燥 对 比 曲 线 
Fig. 4 Comparison of drying curves from one-component 


model, multi-component model and experiment in 70°C 


进入 干燥 未 期 ， 随 着 玉米 籽粒 含水 率 减 小 ， 模 


拟 值 与 实验 值 的 差 值 逐 渐变 大 ， 此 时 模拟 含水 率 变 
化 慢 ， 而 实测 含水 率 变化 快 ， 原 因 可 能 在 于 : 干燥 


实验 末期 ， 在 高 温 下 持续 干燥 较 长 时 间 的 玉米 籽粒 
开裂 ， 玉 米 籽粒 失去 了 籽粒 皮层 阻力 ， 使 籽粒 内 部 
直接 接触 热风 ， 加 快 了 水 分 传递 速率 ， 同 时 玉米 籽 
粒 收 缩 变形 ,水 分 在 玉米 籽粒 内 部 的 传递 距离 变 短 ， 
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相应 的 单位 时 间 水 分 传递 速率 就 变 快 ， 而 在 模拟 实 
验 时 ， 模 型 并 未 考虑 玉米 在 干燥 过 程 中 的 开裂 以 及 
收缩 变形 ， 因 此 造成 了 干燥 末期 含水 率 模拟 值 与 实 
验 值 的 差 值 较 大 的 现象 。 


0.4 、 
实验 什 
a O35 。 多 组 分 
MM 。 单 组 分 
这 一 一 干 基 含 水 率 0.14 
论 0.25 上 
关 
| 
TT 0.15 
0.1 上 1 1 | | 
0 20 40 60 80 100 
干燥 时 间 /min 


图 5 80°C 下 单 组 分 与 多 组 分 模拟 值 与 实验 值 干燥 对 比 曲 线 
Fig. $ Comparison of drying curves from one-component 


model, multi-component model and experiment in 80°C 


对 比 图 4 和 图 5 可 以 发 现 ， 热 风 温 度 越 高 ， 玉 


米 籽 粒 达 到 安全 含水 率 所 需 的 时 间 越 短 ， 干 燥 耗 时 
越 短 。 同 时 ， 可 以 观察 到 热风 温度 80C 时 ， 单 组 分 


模型 与 多 组 分 模型 的 模拟 值 与 实验 值 的 差异 都 较 热 
风 温 度 70'C 时 为 小 ， 这 可 能 是 因为 : 温度 越 高 ， 扩 
散 系 数 越 大 ， 同 时 各 组 分 的 扩散 系数 差异 减 小 ， 接 
近 实 际 值 ， 进 而 模拟 值 与 实验 值 的 差异 变 小 ， 但 考 
虑 到 减少 的 差 值 与 提升 热风 干燥 温度 所 需要 的 能 
量 ， 以 及 高 温 易 降低 玉米 品质 ， 并 不 值得 用 过 高 的 
温度 进行 干燥 ， 以 达到 节能 的 目的 。 同 时 ， 即 使 在 
热风 温度 为 70C 的 时 候 , 多 组 分 的 模拟 值 与 实验 值 
的 差异 也 很 小 ， 说 明 多 组 分 玉米 干燥 模型 较 单 组 分 
玉米 干燥 模型 具有 更 高 的 精度 ， 证 明了 模型 的 可 靠 
生 。 不 同 干燥 温度 条 件 下 ， 玉 米 籽粒 含水 率 子 随 干 
燥 时 间 t(min) 的 变化 规律 如 表 3 所 示 。 

表 3 玉米 干燥 曲线 回归 模型 


Table 3 Formulas of corn drying curves under different 


Pa i 


temperatures 


温度 /C ”模型 关系 式 R? 
Ee a rt 
实验 值 ”X=0.000027 - 0.0044t+0.3905 0.9990 


70 ”多 组 分 了 =0.000027 -0.00441+0.3807 0.9976 
单 组 分 Y=0.00003f2-0.00441+0.3608 0.9805 
实验 值钱 = 0.000037 -0.00521+0.3724 ”0.9963 

80 ”多 组 分 于 =0.00003?-0.00511+0.3657 0.9947 
单 组 分 “和 = 0.00003 已 - 0.0050t + 0.3475 0.9764 


此 外 ， 从 多 组 分 玉米 干燥 模型 的 水 分 分 布 云图 
中 也 可 以 看 出 玉米 籽粒 内 部 的 水 分 变化 情况 。 图 6 
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分 别 为 不 同 温度 条 件 下 ( 左 图 为 70°C , 右 图 为 80'C)， 
玉米 籽粒 模拟 干燥 600s、1800s、3600s 时 的 水 分 分 
布 去 图。 由 图 可 知 ， 随 着 干燥 时 间 的 进行 ， 玉 米 籽 
粒 的 平均 合 水 率 逐 渐 降 低 ; 玉米 籽粒 外 部 含水 率 低 
于 内 部 含水 率 ; 在 不 同 干燥 阶段 , 干燥 温度 为 80'C 
时 的 玉米 籽粒 的 含水 率 较 干 燥 温 度 为 70'C 时 的 低 ， 
即 温度 越 高 干燥 越 快 。 


时 间 =600 动向 : 因 变 量 X (1) 时 间 =600 动向 : 因 变 量 X (1) 
全 0.3904 全 0.3885 


! 


0.4 0.4 


W 0.0401 W 0.0391 


时 间 =1800 动 而 ;: 因 变 量 X (D 时 间 =1800 动 而 : 因 变 量 X (D 
生 0.3526 全 0.3318 


0.2 

0.15 

0.1 

A 0.05 


0 2 
W 0.0401 W 0.0391 


时 间 =3600 动 而 ;: 因 变 量 X 中 时 间 =3600 动 而 ;: 因 变量 XD) 
全 0.2889 人 0.2621 


0.4 0.4 


W 0.0401 W 0.0391 


6 不 同 干燥 时 刻 不 同 干燥 温度 的 水 分 分 布 云图 
( 左 70 'C, 右 80 C) 
Fig.6 Moisture distribution nephogram in different drying 
time of corn kernel ( left: 70°C, right: 80°C) 
4 结论 

(1) 本 文 建立 了 玉米 籽粒 的 多 组 分 与 单 组 分 几 
何 模 型 和 干燥 数学 模型 ， 模 拟 研 究 与 注 层 干燥 实验 
研究 结果 表明 ， 基 于 玉米 籽粒 的 成 分 、 物 理 结构 特 
点 构建 的 多 组 分 干燥 模型 与 单 组 分 干燥 模型 均 可 有 
效 反映 注 层 干燥 试验 曲线 的 变化 趋势 ， 但 多 组 分 干 
燥 模 型 的 模拟 精度 优 于 单 组 分 模型 。 

(2) 干 燥 温 度 升 高 ， 模 型 模拟 结果 与 实验 结果 间 
的 差异 降低 , 但 降低 幅度 有 限 ,热风 温度 为 70C 时 ， 
多 组 分 干燥 模型 依然 具有 较 高 精度 ， 具 有 实际 应 用 
价值 。 
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